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ABSTRACT 

Intensity scintillation has been studied for a He-Ne 
(6328A) laser beam propagating over a 4.05 kilometer hori- 
zontal path in the marine boundary layer. The results 
fall in three areas. (1) The covariance of the logarithm 
of the laser beam amplitude increases significantly with 
decreasing detector aperture diameter for apertures less 
than the first Fresnel zone size. This effect appears to 
be dependent on the magnitude of the refractive index 

structure constant, C^. (2) The dependence of on 

- 4/3 

(height) for unstable meteorological conditions was 

verified using the laser beam in determining C^, • (3) Fre- 

quency analysis performed on scintillation data indicates 
that the increase of power in higher frequencies is more 
dependent on the wind velocity perpendicular to the propaga- 
tion path than on detector aperture diameter for apertures 
smaller than the first Fresnel zone size. 
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I. INTRODUCTION 



The structure of the atmosphere is known to modify a 
laser beam propagating through it. The light may be 
absorbed due to atmospheric gases and pollutants, scattered 
by molecules, electrons or particles, or undergo refraction 
due to turbulence caused by wind and temperature and pres- 
sure gradients. Refraction effects can be further divided 
into beam wander, beam spread and beam scintillation. 

Overall, these effects comprise a significant source of 
errors and distortions in laser communications, radars and 
control systems. 

A. PROJECT HISTORY 

In 1971 an interdisciplinary group of professors at the 
Naval Postgraduate School formed an Electro -Optics group for 
the purpose of pursuing advanced research in a wide variety 
of electro-optic areas. This group has members from Physics, 
Meteorology, Oceanography, Mechanical Engineering, Electrical 
Engineering, and Aeronautical Engineering. A team has been 
investigating aspects of laser beam propagation over the 
ocean surface under a Navy research contract. Previous 
theses by physics students have dealt with investigation of 
the modulation transfer function [Ref. 8] and with verifica- 
tion of the Kolmogoroff turbulence model and comparison of 
C^j as determined by optical methods to that determined by 
micro-meteorological methods [Ref. 6]. This thesis 
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undertakes the investigation of laser beam scintillation in 
the context of aperture-averaging effects, height dependence 
and frequency spectra. Other work being done simultaneously 
involves the determination of laser beam wander and beam 
spread [Ref. 9] . The work thus far has been with 6328 
Angstrom He-Ne lasers. 

B. SCINTILLATION 

1. Description 

Scintillation means "twinkling" or "sparkling" and 
is readily observable if one looks at distant astronomical 
bodies. The twinkling is due to changes in the refractive 
index of the atmosphere. 

A simple model of laser beam scintillation can be 
described by considering the propagation path to be made 
up of small parcels of uniform density (and, therefore, 
constant refractive index). A light ray is bent slightly 
by each parcel, depending on the value of the parcel's 
refractive index, n. A slight variation in path between 
adjacent wave fronts causes interference patterns to be 
formed. The effect at a point receiver is a time varying 
intensity signal. 

2. Theory 

The inherent randomness of turbulence leads to a 
statistical modeling of scintillation. The refractive index 
variations along the propagation path modulate the intensity 
of a beam in a multiplicative manner; thus they modulate the 
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logarithm of the intensity (and amplitude) in an additive 
manner. The variation of the logarithm of the amplitude is, 
therefore, the sum of the random perturbations induced along 
the path. As a consequence of the central -limit theorem, 
the variations of the log-amplitude should follow a normal 
distribution. 

The statistics of the turbulent mixing of the 
atmosphere are provided by Kolmogorof f ' s analysis [Ref. 1]. 

His theory of turbulence allows the mean square difference 
of the refractive index, measured at two points separated 
by a distance d, to be written as 

(1) D n (d) = <[n(r) - n(r+d) ] 2 > 

D n (d) = d 2/3 C N 2 . 

D (d) is called the refractive index structure function. 

2 

is the refractive index structure constant and param- 
eterizes the total amount of energy in the turbulence. 

Tatarski [Ref. 2] has formulated the three-dimensional 
spectral density of the turbulence as 

(2) <f> n (K) = 0.033 C N V 11 / 3 exp[-K 2 /K m 2 ] . 

Here K = 2 tt/£ and is the spatial wave number corresponding 
to Z, the coherence length of a uniform turbulent parcel of 
air. is the spatial wave number corresponding to the 

inner scale of turbulence, Z . This length is the lower 
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limit of the inertial sub-range of turbulence. For values 



of H less than SL , viscous dissipation of energy begins to 
dominate. The outer scale of the inertial sub-range is L q 
and corresponds to the range in which energy is added to 
the turbulence spectrum. The Kolmogoroff theory applies 
only to the region 8, <H<L . The value of H has been 
estimated to be on the order of a few millimeters while L 



assuming that any polarization fluctuations can be neglected, 
and that A<<£ q . X is the wavelength of the propagating beam 
and k now is equal to 2 tt/A. Using Rytov's method [Ref. 3] 
and letting E = exp[if>], (3) becomes 



o 



is on the order of meters. 



The scalar wave equation is given by 



(3) 



? ? ? 

V Z E + knE = 0 



(4) 



V 2 \p + (Vi|mVi|>) + k 2 n 2 = 0 . 



Applying a first order perturbation of 



V = 4> 0 + ^ 



and 



n = n, + 1 



to equation (4) leads to 



(5) 



V 2 ^ 1 + 2 V^ 0 *V^ 1 + 2k 2 ni = 0 . 



8 



The solution is 



•> 



(6) ^(r) = 1/E 0 (r) / 2k n-^ (r^) Eq (r ' ) G (r - r ’ ) dV ' 



Working with equation (6), Schmeltzer [Ref. 4] 

obtained the statistical expression for the log-amplitude 
2 

covariance, C^(d) , which Fried [Ref. 5] applies to a 
spherical wave propagating in a random medium with the result 



Z °° 

(7) C 2 (d) = k 2 /8Tr Re{ / ds / <f> (K 1/2 ) J n (K 1/2 ds/z) 

ha 0 0 n u 



X [1 - exp (Ks (z-s)/iKs]dK 



and 



(8) C £A (d = 0) = C Ia = °- 124k 7/6 Z n/ V 



where 

C 2 = < [ In (A) - ln<A> ] 2 > 



In going from (7) to (8) the assumption <<1 

is made allowing the binomial expansion of the exponential 
and the dropping of higher order terms. This is a weak 
restriction at optical wavelengths. 

The log-intensity covariance is then 



C 



2 

£1 



4 C 



2 

IA 



(9) 



which leads to 



( 10 ) 



res tr 

1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 
9. 

10 . 



11 . 



C 



N 



1.42 V^V 11712 C tI 



The following is a summary of the assumptions 
ictions that apply to this result: 

The turbulence energy spectrum can be described 
by the empirical Kolmogoroff model (equation 1) . 
Polarization fluctuations may be neglected. 

Only the turbulent eddies in the inertial sub- 
range contribute significantly to scintillation. 
The wavelength propagated is much less than the 
inner scale of turbulence, H . 

The fluctuations in the refractive index are 
much smaller than the mean value. 

\ l l/l 4 <<1 . 
o 

Turbulence is isotropic and homogeneous. 

The wave propagated is spherical. 

The detector of the intensity changes is a point 
detector located on the beam axis. 

The first order perturbation technique used in 
the Rytov method adequately models the wave 
propagation. 

The beam is propagated over a horizontal path. 



and 
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II. PROBLEM 



A. APERTURE AVERAGING 

The basic Kolmogoroff model depends on the experimental 
detector being a point receiver. The general effect of having 
a finite detection surface is to reduce the degree of inten- 
sity fluctuations sensed, to reduce the frequency of observed 
changes in intensity and to increase the average signal 
strength. 

1 . Theory 

For atmospheric laser propagation, the amplitude 

1/2 

correlation distance is of the order of (AZ) ' , the first 

Fresnel zone size [Ref. 2] . 

The first Fresnel zone is defined on the collection 

aperture as that radius R at which the geometric path from 

the source to R and from R to the detector is one-half a 

wavelength longer than the axial path. If the diameter of 

1/2 

the receiving aperture is greater than (AZ) ' it will 
contain wave front sections with fluctuations of opposite 
sign so that the overall light flux will fluctuate weakly 
compared to that seen by an aperture with a diameter less 
than the first Fresnel zone size. 

a. Wave-Optics Approach 

Let I(x,y) be the light intensity. The total 
light flux through the aperture is 
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( 11 ) 



P = // I (x } y) dxdy . 

A 

The fluctuations in P are defined as P 1 = P - <P> , and 



(12) P' = // I’(x,y) dxdy 

A 

where I' = I - <I> . 

The mean square fluctuations are 

(13) <P' 2 > = //// <I' (x 1 ,y 1 ) I’ (x 2 >y2' )>dx l dx 2 dy l dy 2 * 

A 1 A 2 

Let Wj = <I' (x 1 ,y 1 ) I ' (x 2 ,y 2 )> = Wj(r) 

where r is the distance between the two points. 

Define F(x,y) such that it is zero outside the aperture and 
1 on its surface. Also define 



x = X-l - x 2 



y = y ! - y 2 



Equation (13) may then be written as 



(14) <P' > = // Wj (x ,y) dxdy // F(x^,y^)F(x^-x,y^-y) dx x dy^ 



Now let 



(15) H(x,y) = // F(x 1 ,y 1 )F(x 1 -x,y 1 -y)dx 1 ,dy 1 
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Evaluating H for a circular aperture of radius R 

F ^ x 1 » y 1 ) = 1 £or x l 2 + y l 2 - r2 
F(x-x 1 ,y-y 1 ) = 1 for (x-x^ 2 + (y-y^ 2 < R 2 . 

The first region is a circle of radius R centered at the 
origin. The second is a circle of radius R centered at 
(Xf.yi) so H(x,y) is the area formed by the intersection 
of the two circles. 

(16) H(x,y) = H (r) = 2R 2 [arccos (r/ 2R) -r/ 2R(1 -r 2 / 4R 2 ) 1/2 ] , r<2R 

= 0, r > 2R . 

Substituting equation (16) into (14) , 

2R 

(17) <P' 2 > = 2 it / W (r) II (r) r dr . 

0 



H(r) equals tt/2 for r = 0, (point aperture), and goes to 

zero for r = 2R. Tatarski [Ref. 2] states that for 
1/2 

R<<(AZ) , the function Wj (r) is approximately equal to 

Wj(0) for 0 r <_ 2R. For this case 

(18) <P' 2 > = (ttR 2 ) 2 <I- I > 2 . 

Tatarski also uses the ratio 



G(R) 



_ Q(R) 
' QTOT 
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where 



Q(R) 



<p» 



<P> 



> 

2 



and 



Q ( o ) = 



<i-i > 2 



The ratio G(R) shows by what factor the fluctuations of the 
light flux through an aperture of radius R are less than 
those through a vanishingly small aperture. 

2R 

(19) G (R) = 4/ttR 2 / W T (r)/W T (0) [arccos (r/2R) 

0 1 1 

- r/2R(l-r 2 /4R 2 ) 1/2 ]rdr . 



G(R) equals unity for R equal to zero and monotonically 

decreases as R increases. 

To accurately evaluate (17), the functions Wj (r) 

and Wj(0) must be known. Fried [Ref. 10] has written G(R) 

2 

(his 0 ) in terms of > the covariance of the log-amplitude 



c jtA^ r ) = ' <£> ] [ £ 0 2 ,y 2 ) - <£>]> 



Where £(x,y) is the log-amplitude and is related to l(x,y), 
the intensity, by 



IO,y) = I Q exp [ (2& (x ,y) ] 



As noted above, £(x,y) is Gaussian random variable. Using 

9 

his data to get values for and assuming different 

2 

values for C HA (0) , Fried calculated G(R) for normalized 
aperture diameters D/ (4z/k) 1//2 . 
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In the limit of the normalized aperture diameter, 
he predicts G(R) should vary as the inverse square of the 
aperture diameter. 

(20) G(R) = D 2 x constant. 

Tatarski [Ref. 2], on the other hand, predicts 

(21) G (R) ~ (D/2(AZ) 1/2 )' 7/3 
1/2 

for D >> 2 (AZ) ' and small values of (see equations (8) 

and (9) ) . 

The main restriction on this theory is that it 
applies to infinite plane waves for horizontal propagation 
paths . 

b. Geometrical Optics Approach 

Young, an astronomer, approached the problem 

using geometric optics [Ref. 11] . For an aperture much 

larger than the inner scale of turbulence or the diffraction 
1/2 

scale, (AZ) , geometrical optics is valid. Assuming the 
aperture is small compared to the outer scale of turbulence, 
the Kolmogoroff isotropic turbulence spectrum produces a 
shadow pattern with modulation power density proportional 

3 

to Z K, where K is the spatial wavenumber. The fractional 
modulation power of K remaining after integration over a 
circular aperture of radius R is 

f (RK) = [2J 1 (RK)/RK] 2 
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where f(RK) is the aperture filter function. 

The total modulation power then is 

(22) P ~ // Z 3 K 1/3 f(RK) K dK d<(> 

= 2ttZ 3 R' 2 / K‘ 2/3 J 2 (RK)dK . 

0 1 

Letting U = RK gives 

(23) P ~ ZV ? / 3 / U-2/ 3 J 2 (U)dU = 0.66 Z 3 R" 7//3 . 

0 1 

Inspection of the integrand of equation (22) shows that it 

falls off rapidly for large K; most of the contribution comes 

from a narrow range of K Z 1/R. Observations of stellar 

3 - 7/3 

scintillation confirm P ~ Z R in the large aperture 

range . 

A surprising feature of this result is for 

R = (XZ) 1 ^ 2 

P - z 3 [CXZ3 1/2 ]- 7/3 = V 7/6 z 11/6 

the result for diffraction-limited turbulence (i.e., very 
small apertures) . This is unexpected for a theory formu- 
lated on geometric optics. 

The restriction on this theory is that R must 
be >>(AZ) 1//2 . 

2 . Published Experimental Data 

Fried, Mevers and Keister made laser beam 
scintillation measurements with 21 collection apertures 
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ranging from one millimeter to one meter in diameter [Ref. 7] . 

These measurements were made with a 6328A He-Ne laser over a 

sloping 8 Km path. The terrain was rolling so that the beam 

had an average height above the surface of greater than 15 

meters. Data accumulation periods were one minute long with 

a time separation between data periods of four minutes. 

2 

Figure 1 is a reconstruction of their data. C^. decreased 

with aperture diameter from 0.438 at 1 mm to about 0.378 at 

2 

10 mm aperture, or about 86 percent of (1 mm). From 

10 mm to 100 mm the drop was from 0.382 to 0.220 such that 



'£A 



(100 mm) ~ 0.557 (10 mm). The data did not dispute 

Fried's theory as presented in Ref. 10 and Section A.l.a. 

above. His data should be evaluated in terms of possible 

2 

sources of error in measurement due to non-s imultaneous 

receiver background compensation, and the assumption that 
2 

was the same over the approximately one and one half 
hours during which the data were taken. 

Young's theory [Ref. 11] has been verified in 
Refs. 12 and 13. This, however, applies to large apertures 
and stellar scintillation observations and would not be 
expected to apply to horizontal laser beam propagation with 
receiver diameters less than the first Fresnel zone size. 

Kerr has reported some limited aperture averaging 
experiments [Ref. 14] . The data consisted of comparisons 
of a 32-cm receiver aperture and a 3 mm aperture at 4880 
Angstroms. Kerr points out that the predicted dependence 
of the irradiance variance is D X for the diameter much 
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1/2 

larger than (A.Z) ' and that it was not observed. Different 



authors assign values of 2, 7/3 or 3 to X. The ratio of log- 

amplitude variances was plotted versus strength of turbulence 

as obtained from micrometeorological measurements. The ratio 

showed an increase from an average of 0.025 to about 0.080 

as the turbulence increased with the largest data spread in 

the saturation region. (See Section II-C for a discussion 

of the saturation region.) In Ref. 15, Dunphy and Kerr 

presented typical ratios in comparison of a 32 cm receiver 

and a 6 mm receiver. The ratios were from 0.35 to 0.70 for 

values of C^ 2 of 5.4 X 10 ^ to 3.1 X 10 ^ 2 . The ratios did 

2 

not always vary consistently with and appear to be 

representative of the saturation regime. 

Kerr summed up the problem of aperture averaging by 
saying, "An adequate theoretical description of receiver 
smoothing has not been given, even for nonsaturated condi- 
tions," [Ref. 14]. It should also be noted that al] of 
the above data were collected from propagation paths over 
land. 

B. DEPENDENCE OF SCINTILLATION ON HEIGHT, THEORY 

The relationship between the refractive- index structure 
2 

function, C^ , and micrometeorological parameters is given 



by 



(24) 




(79 X 10" 6 ^)" 2 C m 2 
rj.2' N 
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2 

is called the temperature structure parameter, P is the 

pressure in millibars and T is the temperature in degrees 

Kelvin [Refs. 6 and 16]. As temperature and pressure depend 

on the heights they are measured at, Wyngaard, Izumi and 

Collins [Ref. 16] developed a semi-empirical theory to 

2 

describe the dependence of on height. From this theory, 

2 

C^j , as determined optically at two separate levels, can be 

written as a function of height. 

2 

C,p may also be written as 

(25) C T 2 = 3.2 N b" 1/3 

where b is the rate of molecular dissipation and N is the 

1 - 2 

rate of molecular destruction of <(T-T) >. Surface-layer 

similarity theory states that statistics of the mean and 

turbulent flow fields, when properly made dimensionless, are 

universal functions of only a stability parameter. By choos- 

2 

ing height h, the kinematic surface stress U* , and the 

d 0 

local mean potential temperature gradient ^ to make the 
statistics dimensionless and the Richardson number Ri as 
the stability parameter, the following relations hold: 

(26) F (Ri) = ^ , and 

1 IT 

U * 



(27) F (Ri) = l l 

hU*O0/3h) 
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Combining (26) and (27) with (25), 

(28) C T 2 = h 4/3 (||-) 2 F 3 (Ri) 



where 



F 3 ( R i) 



3.2 F, 

^T73“ ' 
1 



9 0 

The dependence of and Ri on h must be further 

2 

delineated to establish the relation between Cj and h. 

To do this Wyngaard, et al. [Ref. 16] chose the surface 

9 0 

temperature flux Q in place of and the dimensionless 
parameter j- in place of Ri as the stability parameter. 

L is defined by 



(29) 



-u* 3 t 

kgQ 



where k is Von Karmen's constant and is approximately 
equal to 0.4. 

During the day over land, when heat transfer is generally 
upward, Q is positive and ^ is negative. These are unstable 
conditions, i.e., the buoyant air near the surface tends to 
rise. At night over land, Q is generally negative. Over 
the ocean surface Q may be negative, zero or positive 
depending on the relative temperature influence of the 
ocean. When Q = 0 , ^ = 0 ; these are neutral conditions which 
may be approached under very cloudy conditions. 

A 

To clarify equation (28) , first convert ^ and Ri in 
terms of the new variables, h, U* , Q and L. 
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(30) 



hLL 36 



Q 3h 



h 36 f h^ 

T* 3h " g l l L J 



and 

(31) 



Ri S 2 V 



0 h 

where T*= -fj— • The g’s are functions of r only and are 

u * L 

given in Ref. 17. Then, equation (28) becomes, 



(32) 



C T 2 = T* 2 h“ 2/3 g 3 (£) 



The function g^ is given by 



(33) 



d-u , -2/3 



g 3 = 4.9 [1 -7 (£)]*'*' for 0 > £ 



and 



(34) 



g 3 = 4.9 [1 + 2. 75 (£)] for 0 < £• 



For unstable conditions, with 0 >_ jp and - ^ >> 1, (33) 



~ 4 ( In - 2/3 



becomes g, ~ ^ ( _ r) 



and (32) is then 



(35) 



2 ~ 4 



~ J T* ( - L) 



2/3 , -4/3 



For stable conditions, C^, decreases with height more 
- 2/3 

slowly than h ' . For neutral conditions, the dependence 

-2/3 

approaches h . 

When comparing simultaneous measurements at two different 
heights in the surface layer the following relationship holds, 
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( 36 ) 



C t 2 (1) fh(l)T 4/3 




C T 2 ( 2) |_h ( 2) _ 



For optically determined data then, 



(37) 



C iA (1) C N 2(1) _ r h C2Jl 4/3 



c 2 A (2) VU) u h(ij_ 



under unstable conditions. If conditions are neutral 



(38) 



C t 2 (1) C n 2 (1) _ C^ a 2 (1) fh(2)l 2/3 



C T 2 (2) C n 2 (2) C £A 2 (2) [h ( 1 ) j 



The authors of this theory have pointed out certain 
restrictions to it [Ref. 16]. First of all, this is a sur- 
face layer model requiring local isotropy and may not apply 
at greater heights in the atmosphere. Secondly, the relation- 
ships used hold for long-term, or ensemble, averages of 
turbulence properties. Time averages are only approxima- 
tions to ensemble averages. Thirdly, the theory is limited 
to horizontally homogeneous conditions. This means that a 
given portion of surface, with area large compared to that 
of a typical roughness element, looks like any other portion. 
This is to ensure that the flow structure is in equilibrium 
with the surface. 
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C. SATURATION REGION 
1 . Description 



It is generally accepted in the literature that 
equation (7) , 



C 



2 

SLA 



0.124 k 



7/6^11/ 6 



J N 



2 



does not describe the refractive index structure constant 
for values of 

C?. = < (InA - InA) 2 > > 0.8 



[Refs. 2, 11, 14, 15, and 18 - 21). The effect of saturation 

2 

£A 



2 2 
is seen in comparing C 0A determined optically with deter- 



mined using values of as obtained from micrometeorological 

2 

methods. A one to one relation is followed until (optical) 

reaches about 0.8 where it levels off and then decreases for 

2 

increasing values of (thermal) [Refs. 14, 20 and 21]. 

Tatarski has investigated the effect of using the 

method of smooth perturbations in deriving equation (7) 

2 

concluding that C^(z) << 1 was a limiting factor [Ref. 2] . 



Using a dimensional argument, Kerr deduces the 



relation 



(Maximum) =0.6 



independent of wavelength, for the onset of saturation [Ref. 
14]. 
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2 . Theory 



DeWolf lias developed a theoretical description of 
the behavior of [Ref. 20] . In so doing, he points out 

that the Rytov approximation of smooth perturbations is 
really restricted by the physical condition that ray bending 
does not deflect rays as much as one Fresnel zone off the 
axis . 



The region of turbulence where C 



£ A 



to generally obey equation (7) , with the b 
2 

being for C^ A << 1. The saturation region 



< 1 is considered 
est application 
can be defined by 



1 < K ^/ 3 kz 2 C 2 
1 m Kz L N 



<< 



k 



L 



2 




where K is the spatial wave number corresponding to the 
inner scale of turbulence, k is the optical wave number, 
z is the path length and L q is the outer scale of the inertial 
subrange. In this region the irradiance remains log-normal. 
The dependence of the log-amplitude variance <dX > is calcu- 
lated to be 

(39) <dX 2 > = 0.41 (K m 7/3 z 3 C n 2 )" 1/6 . 

2 

In terms of (for a spherical wave) 



(40) 



<dX 2 > = 0.33 (K m 2 |)' 7/36 (C 2 a ) " 1/6 



The derivation of (39) was based on plane wave theory. In 
order to solve equation (39) must be determined experi- 
mentally, which is difficult to do. 
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Other authors have attempted to model this region, 
but, to date, none have sufficiently satisfied experimental 



observations 



3. 




[Refs . 

VS . 



11 , 



13, and 22 



27) . 



Because of the lack of an adequate theory to describe 
the saturation region it is necessary to correlate measure- 
ments of as determined optically, C^, and as deter- 
mined by micrometeorological methods, C^. One of the 
long range goals of this project is to be able to predict, 
on a large scale, the behavior of a laser beam on the basis 
of meteorological measurements. The mathematical relation- 
ships for the unsaturated region over land are given by 
equations (10) and (24) . 
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III. EXPERIMENTAL PROCEDURE 



A. GENERAL EQUIPMENT DESCRIPTION 
1 . Detector Package 

A detection system was designed and built to function 

under adverse environmental conditions and be adaptable to a 

variety of experimental parameters (Figure 2) . 

The detectors used were standard capacity type PIN-10 

silicon Schottky photodiodes manufactured by United Detector 

Technology, Inc. These detectors provide a linear response 

from 10 ^ watts to 10 ^ watts, have an N.E.P. (lKHz - .85p - 
- 1 2 

1Hz BW) of 10 watts, a rise time of 10 nanoseconds and 

a dark current of 0.5 microamperes. The units were operated 

in the photovoltaic mode with zero bias voltage and a large 

load resistance (typically 500,000 ohms). The active surface 
2 

area (1.25 cm ) was found to be non-uniform, somewhat 
reminiscent of photomultiplier tubes. A typical plot of 
relative voltage output versus position of a focused 6328 A 
laser beam is shown in Figure 3. 

A 31 cm focal length lens was used in conjunction 
with a thin plate aperture which limited the field-of-view, 
and a 2.54 cm field lens to focus the input aperture on the 
photodiode. The field lens was necessary to remove a source 
of possible error due to either beam wander effects or 
tracking errors. The lens eliminated any motion of the 
image on the detection surface and thereby ensured a 
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constant responsivity from the detector. The field limiting 
plate aperture was placed at the focal point of the objective 
lens. The aperture was removable so that its size could 
be changed to meet the particular requirements of the experi- 
ment undertaken. 

Placed immediately in front of the objective lens 
was a 10 A filter centered at 6328 A with a 65 percent peak 
transmission coefficient. The narrow band filter, along with 
the field limiting aperture helped reduce the effect of back- 
ground radiation. The filter mounting limited the maximum 
effective receiver aperture to a 23 mm diameter. In front of 
the filter was an 18 leaf iris capable of being adjusted 
from 42mm to 3mm. 

The main feature of the detection system was its 
method of background compensation. From Figure 2 it can be 
seen that the detection package was made up of two separate 
receiving paths. The laser signal was received through one 
path and the background radiation capable of being trans- 
mitted through the narrow band filter was received through 
the other. The field limiting aperture in the background 
sensor was horizontally offset from the optical axis so 
that it would not see the laser signal. A typical offset 
was such that the areas seen by the two receivers were 
separated by 1-1/2 meters at 4 kilometers. The placement 
of these apertures was tested in the laboratory over a 
path length of 280 meters prior to use in an experiment. 

The electronics system was then given simultaneous readings 



27 



of the laser signal and the background radiation from a 
field of view as nearly equivalent to that of the laser 
receiver as possible. 

2 . Electronic Components 

A block diagram of the basic electronic components 
used is depicted in Figure 4. The voltage signals from the 
detector package were sent to a low noise differential pre- 
amplifier which subtracted the background radiation from 
the laser radiation, amplified the signal and sent it through 
a circuit box to the log-amplifier. At the circuit box the 
amplified signal was monitored using a microvoltmeter, and 
a calibration signal could be fed into the system when de- 
sired. A bias circuit was also built into the signal line. 

By chopping the signal at the laser every 60 seconds, for a 
10 second period, the operator could see the amplified dif- 
ference between the two detector signals and set the differ- 
ence equal to zero. It should be noted that during daylight 
hours it was necessary to make adjustments to the bias con- 
trol. On many occasions, the difference in background could 
be observed to drift, that is, the background radiation 
changed over short periods of time. This fact was the 
initial reason for using a simultaneous background compen- 
sation system as opposed to first recording the laser signal 
for a period of time and then recording the background, 
followed by a deconvolution of the two sets of data. 

The log-amplifier was capable of accepting signals 
of 1 millivolt to 10 volts. The output was proportional to 
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the logarithm (base ten) of the signal so that input signals 
of 10 mv, 100 mv and 10 v corresponded to outputs of 200 mv, 
400 mv, and 800 mv, respectively. This output was then sent 
to a wide band power amplifier where it was amplified by a 
factor of ten. From the power amplifier the signal was sent 
to an oscilloscope, where it was monitored, and then on to 
an F.M. conversion instrumentation tape recorder where dt 
was stored. 

The final data processing was done back in the 
laboratory where the tape recorded signal was sent to a 400 
channel pulse height analyzer. Half of the analyzer storage 
was used to store the applicable calibration signal while the 
other half was used to accumulate and display the distribu- 
tion of the scintillation signal. The sampling frequency 
was set to 16 KHz by an external trigger pulse generator. 

The laser used for the aperture averaging experiment 
and for the lower level of the two-height experiment was a 
He-Ne Metrologic Model 920 providing an effective output 
power of 2.3 milliwatts at 6328 A. The laser was found to 
be stable after a half hour warm-up period. A 6mm focal 
length microscope lens system mounted on the laser output 
aperture imaged the beam on an 4cm focal length converging 
lens. This converging lens was adjustable in horizontal 
position allowing the beam to be diverged or focused at 
will. This provided an effective transmitter aperture of 
6 mm. The laser with its associated lens systems was then 
mounted on a four inch channel. This channel, in turn. 
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was mounted on a base constructed of 3/4 inch aluminum 
plate in such a manner that an operator could make fine 
adjustments in train and elevation. 

An additional laser was necessary for the two- 
height experiment. This laser was He-Ne, 6328 A, model 
S-100 manufactured by C-W Radiation and provided 0.8 milli- 
watts power. It was fitted with an adjustable beam spreader 
and was mounted on a heavy duty camera mount. 

B. APERTURE AVERAGING 
1 . Propagation Path 

The propagation path was 4.05 kilometers in length, 

98 percent of which was over the ocean surface at a height of 
4.9 meters (Figure 5). The laser was located on the ground 
level of the Monterey Holiday Inn about three meters inland 
from a 2.5 meter high vertical sea wall. The beach sloped 
gently from the sea wall about 40 meters to the ocean. The 
receivers were located in close proximity to a three-story 
building at Hopkins Marine Station. Here, the land sloped 
gently up from the ocean about 30 meters to the receiver 
site. 

Data were taken on two days, 22 October and 
27 October, 1973. The propagation path generally was perpen- 
dicular to the ocean wave flow. On these two days the average 
wave height was about one foot with three foot swells. Both 
days were sunny with no cloud influence. The first day had 
strong winds (20 knots) from the Southwest while the second 
day was calm with variable light winds (0-5 knots) . 
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2. Data Accumulation 



An earlier trial experiment across the same path 
indicated the need for multiple aperture data. The four 
apertures used did not provide a readily defineable relation- 
ship so it was decided to augment the small end of the 
selection with two additional apertures. Pieces of aluminum 
stock were machined to the outside diameters of 15, 9, 6.3, 
4.8 and 3 mm with the sixth aperture being the 23 mm effec- 
tive maximum aperture of the narrow band filters. All 
apertures were less than the first Fresnel zone, 5.08 cm. 

To change from one size aperture to another, the desired size 
aperture gauge was placed in the opening of the irises, 
the irises closed on them and the gauge removed. The move- 
ment of the irises was tight enough to ensure that there 
was no change in the opening once the size was set. 

The major problem discovered in the trial experiment 
was that separate aperture data accumulation periods could 
not be correlated without the possibility of significant 
errors due to the change in C^ T . To circumvent this problem 
another detection package was required to provide a continu- 
ous measurement of C^, at a fixed aperture. The detector 
package described in Ref. 6 was modified to match the 
maximum aperture size of the variable aperture unit and a 
suitable field limiting aperture plate was machined to match 

the other unit. The two packages were tested in the corri- 

2 

dor laboratory to see if they would yield identical 
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using the 23 mm apertures. The systems agreed to ±1.3 
percent under weak turbulence conditions by placing them 
side by side. 

The data taken on 22 and 27 October, then, consisted 
of a continuous monitoring of with a 23 mm aperture while 
varying the aperture size on the other detector package. Each 
complete data sequence started and ended with a run at 23 mm 
on the variable aperture package, taking about 14 minutes. 

Each run was of two minute duration with the laser signal 
being chopped every minute for background bias adjustment. 

The signals from the receivers ran through parallel electronic 
components to the E.M. tape recorder where they were simulta- 
neously taped on different channels. An additional voice 
channel was utilized to make pertinent observations, in- 
cluding time calibrations for each run. Six data sequences 
were taken on 22 October and seven on the 27th of October. 

3 . Data Reduction 

Back in the laboratory, the taped signals were fed 
into an oscilloscope for monitoring purposes and from there 
into the pulse height analyzer, PIP 400. The first signals 
recorded during each experiment were the calibration signals. 
These signals were 50, 100 and 200 millivolts at the input 
side of the log-amplifiers. They were stored on 200 channels 
at the PIP 400 in such a way that they could be superimposed 
over the log-intensity distributions after each run was 
accumulated on the other 200 channels. 
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The intensity distributions were log normal, as 

predicted by theory [Ref. 28]. The desired quantity was 

2 2 

, the variance of the [In (I/<I>)] distribution. The 

quantity displayed on the PIP 400 was 

(41) alog^pl = 0.434 alnl = blnl 

where a and b are constants. The full width at half maximum 
of a normal distribution is X ££ / 2 = 2.345 a. Therefore, 

(42) a = 0.425 X £1/2 . 

Both a and X ££ ^ 2 are measured in units of whatever 
is distributed, in this application blnl. Now 

(43) <(bln ^) 2 > = b2< Cln <|>) 2 > = b 2 C^ . 



To get a in the right units, a scaling factor, S, based on 
the known calibration signals was used. 



(44) 



b [lnl 2 - lnl £ 1 
X 2 -X £ (cm) - 



Then the variance of the distribution could be found by 
measuring X ££ ^ 2 in cm and using 



(45) 



a = 



0.425 bln(I 2 /I 1 )X £1/2 (cm) 



X, 



X £ (cm) 



which in turn yields 
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and finally, 



b [0.425 lnCJ 2 /I 1 )X fl/2 ] 

lx 2 - i' 1 'J 



b C 



n ’ 



(46) 



'0.425 ln(I 2 /I 1 )X fl/2 (cm) 



,2 



'LI 



X, 



XT (cm) 



Each distribution, with the calibration signals 
superimposed, was photographed. was measured off 

the photograph using a vernier calipers. Maximum error 
introduced by this technique was ±2.7 percent. was 

computed using equations (10) and (46) . 

C. TWO LEVEL EXPERIMENT 

1 . Combined Meteorological and Optical Effort 

The purpose of this experiment was to test the theory 
presented in Section II.B. over the ocean surface. The same 
propagation path as described above was utilized for the 
lower level. The" 2. 3 mw laser and its corresponding re- 
ceiver were 4.9 meters above the ocean surface. A second 
system using the 0.8 mw laser measured scintillation at a 
height of 12.3 meters, directly above the lower path. 

The Research Vessel ACANIA was fitted with the 
meteorological instrumentation and anchored about 15 meters 
South of the mid-point of the propagation paths. Data were 
taken from 1943 to 2221 (PDST) on 18 October. The wind 
was from the South at 3 - 5 knots ; waves were about two 
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feet high with two foot swells. For most of the experiment 
the area was covered with low clouds. 

2 . Data Accumulation and Reduction 
a. Meteorological 

This portion of the experiment was performed by 
members of the Meteorological Group. Meteorological instru- 
ments and measurement techniques were essentially the same 
as those reported in Ref. 29. An instrumentation mast was 
installed on the bow of R/V ACANIA. Two amplifier bridge 
circuits, (Sylvania GTE Thermosondes), each with two 60 
ohm platinum wire sensors were mounted at 4.525 meters and 
8.41 meters above the water line. could then be obtained 

by measuring changes in the resistance difference between 
the two wire sensors spaced 10 cm apart. The resistance 
difference is proportional to the temperature difference. 

A voltage signal proportional to AT (rms) was recorded on a 
strip chart as a function of time. 

In processing the data, raw values of AT (rms) 
were abstracted from the recorded chart. These values were 
then corrected for calibration errors in the nominal resis- 
tance of the platinum wires. 

C T was then calculated using 

C^, = true AT (rms)/r'*'^ 

where r is the probe separation of 10 cm. From this, 
was found using equation (24) . 
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Stability conditions were estimated from observa- 
tions of the gradient formed by the sea water temperature 
and air temperature sampled continuously throughout the 
experiment. Quartz thermometers were placed at the lower 
level and on the ocean surface, 
b. Optical 

The equipment used to receive and process the 
data was the same as for aperture averaging with one addi- 
tion. All the electronic components were located at the 
lower receiver site to limit the number of operators to 
one. This required the use of a long coaxial cable from 
the upper receiver to its differential amplifier. To ensure 
no losses were encountered, a cable driving cathode- follower 
circuit was built using field effect transistors. The 
circuit performed the function of subtracting the background 
signal from the laser signal as well as drying the long 
cable . 

The field limiting aperture plates were machined 
to cover a 9.15 meter diameter circle at the 4.05 Km distance. 
This, together with 3.7 meter beam diameter sizes at the 
receiver ensured that each receiver was only seeing the out- 
put of its respective laser. 

Data were accumulated for an average of three 
minutes per run. All data reduction was accomplished in 
the laboratory using the same techniques as described in 
Section III.B.3. 
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D. FREQUENCY SPECTRUM ANALYSIS 



Representative samples of data were chosen to compare 
the relative power spectra of scintillation. These samples 
include upper and lower level data from the two-height 
experiment, comparison of the two detector systems when side 
by side (12 cm apart) , and a complete aperture averaging 
sequence . 

The recorded scintillation signal was sent through a 
wave analyzer with a band-width of 10 Hz. The lower 
frequency limit of the analyzer was 20 Hz. Each signal was 
filtered at 10 Hz intervals from 20 to 100 Hz, at 100 Hz 
intervals from 100 to 1000 Hz, and at 1500, 2000, 3000, 

4000 and 5000 Hz. The output was sent to an X-Y recorder. 
The power was determined by integrating the output curves 
with a planimeter. 
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IV. EXPERIMENTAL RESULTS 



A. APERTURE AVERAGING 

1 . Preliminary Investigation 

Initial investigation of the effect of aperture 
2 

averaging on using diameters smaller than the first 

Fresnel zone was conducted in the laboratory. The C-W 
Radiation 0.8 mw 6328 A laser beam, slightly diverged, was 
sent down a long corridor and reflected back to the receiver 
by a plane mirror. Total path length was 286 m. Turbulence 
was provided by seven large forced draft ventilators which, 
admittedly, was not representative of turbulent flow in 
the open atmosphere, but nevertheless, did produce a log- 
intensity distribution that proved to be normal. 

Each sequence of measurements contained data taken 
using apertures of 3.5, 9, 18 and 23 mm. Relative values of 
were calculated by using obtained with the 3.5 mm 

aperture as the normalizing factor. The results of five 
data sequences and a composite plot are shown in Figure 6. 
Each data point is an average value taken from five or more 
30 second sampling periods. The degree of aperture averaging 
found here indicated that further investigation on a full 
scale path was necessary. 

2 . Monterey Bay Experiments 
a. 22 October 1973 

As previously mentioned, this day was somewhat 
unusual in that the winds were quite strong and were from the 
South, i.e., over land. 
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The control receiver system, operated with a 
fixed 23 mm aperture, was used to continuously record data 
from 1555 to 1751 (PDST) . Figure 7 is a plot of the values 
of obtained when the variable aperture system was set 

to 23 mm vs. time. Only one data run provided identical 
values of from both, receiver systems. For the remainder 

of the runs, values obtained with the variable aperture 

system were as much as 12 percent lower than the control 
system values. The average deviation was on the order of 
6 percent lower. Some factors that may be responsible for 
this difference are beam wander, instrumental error and 
position of the receivers with respect to the beam axis. 

Assuming that the calibration set in the 
laboratory was still valid leads to an investigation of 
a combination of beam wander and position relative to the 
beam axis. (Beam wander was measured at approximately 15 
minute intervals by Beall and is reported in Ref. 9.) 

The intensity distribution of the beam in the absence of 
any turbulence can be considered to be Gaussian. If a 
detector is placed in the beam and the beam is allowed to 
wander randomly, the intensity differences sensed will 
depend on the magnitude of the beam wander. If the detector 
is on the original axis, the highest intensity sensed will 
be when the beam center is coincident with the detector and 
the lowest intensity will be when the beam is at its maximum 
excursion from the original axis. If, however, the detector 
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is placed at a point further off the original axis than the 
maximum beam wander, the maximum intensity sensed will be 
at one limit of the maximum beam wander and the minimum will 
be at the other. The intensity difference, then, will be 
much larger than that sensed when the detector is on the 
beam axis. When propagating in a turbulent atmosphere, 
these intensity differences are coupled with those caused 
by the medium inhomogeneities and increase the apparent 
scintillation. 

On this day the position of the receivers in 
the beam was found to be about 0.35 m off center. The beam 
diameter was 2.74 m and the distance between sensors was 
0.12 m. The control (i.e., reference) system was further 
away from the beam axis and, as shown in Figure 7, always 
had the larger value of . 

Ho [Ref. 30] has shown theoretically that the 
variance of the log-amplitude fluctuations, C^, increases 
with radial distance from the axis of a laser beam when the 
aperture is located in the near field. The 4.05 km propaga- 
tion path and 3 mm radius transmitter aperture placed the 
receivers in the far field. A similar effect to that 
described by Ho was, however, observed in the corridor labor- 
atory under far field conditions. This effect, amplified 
by beam wander, appears to explain qualitatively the dif- 
ference in C f A values obtained. 

The purpose of the control system was to 
continuously monitor the changes in the atmosphere so that 
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the values of obtained with the smaller apertures could 

be compared to a known reference. Because identical values 
of C ^ were not obtained with the 23 mm apertures at the 
beginning and end of each sequence, a linear interpolation 
was carried out to obtain baseline data at intermediate 
times. This interpolation process consisted of comparing 
the difference in values observed at the beginning of 

a sequence to the difference in values observed at the 

end of the sequence. Define these values as AC^(t-^) and 
^ZA^Z^ ’ an ^ F as t * le rate change of AC^(t). Then 



F = 



AC ^A ( ' t l' ) " AC £,A ( - t 2' ) 



The interpolated value of C^^(23mm) for the variable aperture 
system at a time t in the interval t^ <_ t _< would be 
equal to 

C £A^ + f AC £A^ t l' ) " F( ' t-t l^ 



where C^(t) is the value of obtained with the control 

system at time t. 

Figures 8 and 9 depict the relationships between 
C £A CD) as seen by the respective apertures and C^(23 mm) as 
determined by the interpolation process. These figures 
show that C^(D) does not consistently increase and decrease 
with C^(23mm) and that the value of C 0A (D) - C 0A (23mm) is 



l A 



J ZA 



rather consistent for a given aperture. 
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Figure 10 is a plot of the differences between 
the values of C^(D) for the respective apertures and 
as interpolated from the 23 mm aperture measurements. The 
curve rises with decreasing aperture size until 4.8 mm is 
reached. The average value for the 3 mm aperture was slightly 
less than for the 4.8 mm aperture, 
b. 27 October 1973 

This day was sunny with very light and variable 
winds. The same procedures were followed as above. Figure 31 
is a plot of the 23 mm values of as determined by both 

systems. Once again, the two systems did not yield identical 
values. The control receiver system was closer to the average 
position of the beam axis. The difference in the plots may be 
qualitative] y explained using the same arguments as in the 
previous section and adding the fact that the overall value 
of was larger. (Beam wander measurements were not taken 
on this day.) 

Some interesting visual observations of beam 
wander were made. The laser was first sighted in on the 
receiver array around 1200 (PDST) but kept drifting vertically 
downward. Because of the laser elevation adjustment mechanism, 
any vertical drift due to mechanical slippage would have 
caused the beam to move upward; therefore, the beam motion 
was due to the atmosphere. The rate of drop was about one 
meter per 10 minutes. The beam stabilized from about 1330 
to 1530 (PDST) at which time it started to drift upward at 
about 0.3 meters per 10 minutes. A possible explanation 
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would be large scale density changes, or more likely density 
gradients. The sunny calm conditions were conducive to the 
formation of a non-dissipating warm air mass over the cooler 
ocean surface. This could form a density gradient with the 
density greater nearer the surface. The beam would then 
tend to bend downward. As the air above the surface cooled 
down later in the day, the beam would move in the opposite 
direction. 



Figures 12 and 13 show the relationship between 
the aperture averaged system and the extrapolated 23 mm 
values. Figure 14 is a plot of the difference between the 
extrapolated 23 mm values and the values obtained using the 
various apertures versus aperture size. The shape of the 
curve is very similar to the one obtained from the data of 
the 22nd (Figure 10) . The peak of the curve in Figure 14 
is 2.6 times as large as the curve in Figure 10 indicating 
that the aperture averaging effect for apertures smaller 
than the first Fresnel zone may be functionally dependent 



c. Composite Data 

Figures 10 and 14 show that the largest relative 
values of were obtained when the 4.8 mm aperture was 
used. Based on this, it was assumed that the 4.8 mm aperture 
approximated a point receiver better than the other apertures 
used. 



In order to find a relationship between aperture 
averaging and C^, an interpolation process, similar to that 
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carried out on the 23 mm apertures, was undertaken to 
establish a continuous set of values for as seen by the 

4.8 mm aperture. The ratios of (D) (4 . 8 mm) were 
calculated, averaged for each collection aperture each day, 
and plotted as a function of aperture diameter in Figure 15. 
Values of C^^(D) /C^ (1mm) obtained from curves published by 
Fried, et al. [Ref. 7] and reproduced as Figure 1 are 
plotted for comparison. Linear regression lines were calcu- 
lated and plotted for all three data sets. The slope of the 
regression line for the 27 October data was steeper than 
those corresponding to 22 October and Fried's data. The 
lines cross in the vicinity of the point C^^(D)/C^^ equal 
to 0.93 and aperture diameter equal to 10.6 mm. Figure 16 
is a plot of the absolute values of the slope ^ vs as 

seen by the 23 mm apertures for the regression lines plotted 
in Figure 15. The line was calculated using the least square 
method. (Fried's data were adjusted to apply to the 4.05 km 
path.) Using this limited amount of data a point collection 
aperture value of could be estimated in the following 

manner. Entering Figure 16 with the observed value of 
C^(23 mm) yields the inverse slope to be used in Figure 15. 
Applying this slope in such a way as to cross through the 
empirically observed point of intersection of the regression 
lines will give a value for the ratio C^(23 mm)/C^^, which, 



in turn yields C 



a A' 
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It is obvious that additional data are needed 
to formulate a reliable mathematical relationship between 
aperture diameter and C ^ A or C^. 

d. Frequency Spectrum 

Clifford [Ref. 31] has calculated the theoretical 

temporal frequency spectrum, W(f), for a spherical wave. He 

solved the plane and spherical wave asymptotic forms in 
2 

terms of C 0A , f and ft 
£A’ o 

where f Q = v/ (2t AL) 3 , 

v is the wind velocity perpendicular to the propaga- 
tion path, and 

ft = f/f . 

' o 

For a plane wave 



(47) 



W 



(f) = 0.44 (C £A )^ (^-)(l + 0.27ft 4/3 ) } ft « 1 



(48) 



W(f) = 1*14 (C^ A ) p (^-)^' 8/3 ,^ >> 1 



In the spherical wave case 



(49) 



W(f) - 0.245 (C £a ) 3 (j-) [l + 0.119ft 4/3 ], ft << 1 

o 



(50) 



W(£) = 2.81 (C £A )^ (^-)ft' 8/3 

o 



ft >> 1 



These solutions are valid for a vanishingly small aperture. 

Tatarski [Ref. 2] calculated the effect of 
receiver aperture averaging on the frequency spectrum for a 
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plane wave model. This effect shows up in an additional 
term multiplied on the kernal of the integral, which corre- 
sponds to a partial suppression of the normalized fluctua- 
tions. This factor is 



where f is the frequency, v is the wind velocity perpendicular 
to the propagation path and K is the spatial wave number. For 
f >> v/2rrR, this term is small for all K. For high frequen- 
cies the aperture averaging weakens the effect of inhomogene- 
ities on all scales. If f << v/2irR, the factor is close to 
unity for K << ~ and becomes small for K >> jy . The contri- 
bution of large-scale inhomogene j ties in the low-frequency 
part of the spectrum is virtually unaffected, whereas all 
scales which are small compared to the aperture radius are 
suppressed . 

Relative power values for each aperture were 
obtained by dividing the power found at a particular fre- 
quency by that found at 20 Hz. From equations (49) and (50) 
it can be seen that the ratio W(f^)/W(f 2 ) for a vanishingly 
small aperture is a function of f q , f ^ , and f 2 only. The 
aperture dependence is contained as a multiplicative term 
in a first order Bessel function as shown in (51). The 
dependence should not be lost in taking the ratio 
lV(f 1 )/W(f 2 ) ; therefore, when comparing relative power vs. 



2 J 



1 




(51) 



r/k 2 + 4Tr Z f Vv 2 ' 
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frequency plots, a difference should be detectable due to an 
aperture change. If the data samples were not taken simulta- 
neously, changes in v could also contribute to any differences 
observed. Estimates were made of v because the propagation 
path could not be properly instrumented. 

On 22 October, v was estimated to be 16 knots 
(8.24 meters/second). Defining 

(52) f c = v/2ttR , 

then for this day f = 1.31 (1/R) . Here 1^(6. 3 mm) = 417 Hz 
and f (23 mm) = 114 Hz. Figure 17 is a comparison plot of 
relative power vs. frequency for data obtained with a 23 mm 
aperture and for data simultaneously obtained with a 6.3 mm 
aperture . 

For f > 500 Hz, the relative power curves in 
Figure 17 are converging. For most frequencies sampled, the 
23 mm aperture receiver showed a smaller contribution to 
relative power. This is in agreement with the aperture 
dependence and with the previous statement that all scales 
of inhomogeneities which are small compared to the aperture 
radius are suppressed. 

On 27 October the winds were light and variable. 

If an average value of 1 knot is assumed for v, f = 0.082 
(1/R) and f Q = 4 Hz. The value of f ranges from 53 Hz 
for a 3 mm aperture to 7 Hz for a 23 mm aperture. Figures 
18 and 19 are plots of the relative power vs. frequency 
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for the sequence of apertures. The plots are similar in 
shape to Figure 17. 

Figure 18 was plotted from data received by 23 mm, 
9 mm and 4.8 mm apertures. The values of f for v equal to 
1 knot are 7 Hz, 18 Hz and 34 Hz respectively. For f < 90 Hz, 
variation with aperture size of the relative power in any 
frequency band is unpredictable and is probably due to slight 
variations in v for the different sample periods. Sample 
periods for the different apertures were four minutes apart. 
For f > 90 Hz, the smaller apertures show a larger contribu- 
tion to relative power and the values tend to converge on each 
other as they approach zero at about 500 Hz. This is in 
agreement with Tatarski's model for f >> f . The temporal 
variations in v would have less effect in this region. 

Figure 19 data were received by 15 mm, 6.3 mm 
and 3 mm apertures. The corresponding values of f for v 
equal to 1 knot are 11 Hz, 26 Hz and 55 Hz. The same be- 
havior of the relative power vs. frequency is observed as 
in Figure 18 except that the region where ordering of 
W(f)/W(20Hz) is seen appears to be f > 300 Hz. It would be 
expected that this region would begin at a higher frequency 
due to the requirement that f >> f c and the values of f 
for this graph are generally larger than for the apertures 
used in Figure 18. 

Figure 20 is a comparison of relative power vs. 
frequency as received by the 6.3 mm aperture on 22 October 
and 27 October. The values corresponding to 22 October are 
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seen to be consistently larger than those of 27 October. 

This difference is due to the significantly larger value of 
v for 22 October. A comparison of this figure with Figures 
18 and 19 shows that the effect of v on relative power is 
greater than that of aperture size within the limits of 
these experiments. 

Figure 21 is a plot of relative power vs. 
frequency as seen by the 23 mm apertures of the variable 
aperture system and the fixed aperture control system. The 
control system shows slightly lower values of relative power 
for most of the frequencies sampled. 

B. DEPENDENCE OF SCINTILLATION ON HEIGHT 
1 . Optical Measurements 

Optical measurements were taken from 1944 to 2219 
on 18 October 1973 (PDST) . The plot of C^ vs. time for the 
two levels is shown in Figure 22. Both receivers utilized 
a 23 mm aperture. The values obtained at the upper level for 
the first 40 minutes are not considered to be representative 
because of difficulties in maintaining the beam axis on the 
receiver. The large gap in upper level data from 2023 to 
2105 was due to making beam adjustments. The data from 2105 
to 2206 shows a definite relationship between the two levels. 

Table 1 displays the average values of C^ for three 
different averaging periods and their relationship to the 
heights of the two systems. Period (a) was from 2105 to 
2137; (b) from 2140 to 2206; and (c) from 2209 through 
2221 (PDST). 
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TABLE 1 





Period (a) 


Period (b) 


Period (c) 


iA at h (1) = 16 ft. 


0.360 


0.349 


0.271 


£A at h (2) = 40 ft. 


0.192 


0.192 


0.204 


(h(2)/h(l))' l/3 


3.44 


3.44 




(h(2)/h(l)) 2/3 






1.85 


La)/ C L (2) 


3.51 


3.30 


1.7 7 



As discussed in Section II. B. , during stable 

2 2 

meteorological conditions C^, (and therefore ) decreases 

- 2/3 

with height more slowly than h . Under neutral conditions 

- 2/3 

the dependence approaches h and under unstable conditions 

varies as h . For stable conditions with data taken 
simultaneously at h(l) and h(2) 



(38) 



C T 2 W m C N 2 ^^ = C IaM 

C t 2 (2) ' C N Z (2) C £ A ( 2 ) 



h(2)l 



2/3 



h U) J 



For unstable conditions 



(37) 



c Lw = Vd) 



h ( 2 ) 



h (1) 



J 



4/3 



For the results displayed in Table 1, period (a) 
shows that (h ( 2) /h (1) ) agrees with C^ A (1) /CJ A (2) within 
2 percent; for period (b) the agreement is 4 percent. For 
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period (c) (h(2)/h(l)) 2 / 3 agre es with C^ A (1) /C^ A (2) to 

within 4 percent, as would be expected for neutral 
conditions . 

Table 2 shows the same data with aperture averaging 
corrections added using Figures 15 and 16. 



C M at 
C £A at h (2) 

(h(2)/h(l)) 4/3 

(h(2)/h(lj) 2/3 

c\ k a)/ c \m 



TABLE 2 

Period (a) 
0.442 
0.226 
3.44 

3.81 



Period (b) 
0.425 
0.226 
3.44 

3.53 



Period (c) 
0.324 
0.240 

1.85 

1.82 



Period (a) shows agreement to 9.6' percent with the predictions 

for unstable conditions. Period (b) is within 2.5 percent of 

agreement with unstable conditions and period (c) is within 

1.5 percent agreement with neutral stability conditions. 

2. Optical Frequency Spectra 

Figure 23 is a plot of relative power vs. frequency 

as seen by the two levels at 2134 (PDST) 18 October. The 

plot is very much like that of Figure 21 where the two 

detection systems were placed side by side. From equations 

(49) , (50) and (51) , the theoretical temporal frequency 

2 

spectrum, W(f) is a function of C^ A , f , f and R. The 
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ratio W(f)/W(20Hz) varies only with f and f. For a given 
frequency, the relative power values are dependent on the 
effective magnitude of v at the sensor level. Figure 23 
indicates that the value of v must have been comparable at 
the two levels on 18 October and less than v on 22 October. 

3 . Meteorological Measurements 

As defined in Section II. B., when the surface 
temperature flux, Q, is positive, the air mass is unstable. 

Q equal to zero corresponds to neutral conditions and Q 
negative corresponds to stable conditions. The temperature 
difference between the sea water and the air at the lower 
level, 4.525 meters above the surface, is plotted in Figure 24. 
Conditions were unstable from 2105 to 2200 when the ship left 
the area. Although quantitative measurements of Q after 
2200 were not available, low lying coastal clouds were 
setting in throughout the evening with a ceiling of 200-300 
feet by 2220. The increasing cloudiness would tend to 
change conditions toward neutrality and the corresponding 
-2/3 relationship. 

A malfunction in the lower level temperature probe 
system prevented the comparison of the two level optical data 
with meteorological data. The difference of 3.9 meters in 
the upper levels of the optical and meteorological sensors 
prevented a point by point comparison of values. 
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V. SUMMARY AND CONCLUSIONS 



A. APERTURE AVERAGING 

Significant aperture averaging was observed in an 
enclosed laboratory for apertures less than the first Fresnel 
zone size. Two field experiments, conducted under different 
turbulence conditions, indicated aperture averaging, for small 
apertures, is dependent on C^. Data collected by Fried, 

Me vers and Keister [Ref. 7] lends support to this conclusion. 
Both field experiments show that a 4.8 mm aperture can 
approximate a vanishingly small receiver and should be used 
whenever possible. None of the existing theories adequately 
describe the observed data. 

Application of a primitive method of predicting C^(0) 
from C ^ observed by a 23 mm collection aperture appeared 
to be successful for the two level experiment data. Addi- 
tional aperture averaging data is required to make accurate 
predictions of C^,^(0). The capability to make such predic- 
tions would greatly reduce the intensity per unit area 
required to collect good scintillation data. 

The observed frequency spectrum is in general agreement 
with the predicted trends in Clifford's model for a spherical 
wave [Ref. 31] and Tatarski's aperture averaging modifica- 
tion term as applied to a plane wave [Ref. 2] . The effect 
of smaller apertures seeing more relative power in the 
higher frequencies than larger apertures was observed. For 
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the range of apertures used (3mm to 23 mm) the effect of 
aperture size on relative power at a given frequency was 
seen to be much less than the v dependence. An increase in 
v shows up as an increase in relative power in the higher 
frequencies . 

B. DEPENDENCE OF SCINTILLATION ON HEIGHT 

2 

Wyngaard, Izumi and Collins [Ref. 16] predicted C^ 

was dependent on h f or unstable meteorological conditions 

-2/3 

and on h for neutral meteorological conditions in the 

context of long time averages, near the surface, over land. 
The data presented in Tables 1 and 2, though sparse, are 
in excellent agreement with these predictions and indicate 
that the semi-empirical theory is applicable to propagation 
paths over the ocean surface. 

C. EXPERIMENTAL PROCEDURE 

Accumulation of scintillation data for the aperture 
averaging experiments and the two level experiment pointed 
out some potential sources of error. Beam wander makes a 
significant contribution to the scintillation observed by a 
detector. If meteorological measurements are to be used to 
predict the amount of scintillation a laser beam experiences, 
an additional term should be found to compensate for the 
effect of beam wander. 

The value of C^, is strongly dependent on the height 
of the propagation path. This fact can be used to obtain 
large ranges of C^ on any given day. When making comparisons 
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of C^q and C>,^, care should be taken to ensure all optical 
components and meteorological sensors are on the same level 
because the - 4/3 and - 2/3 scaling lav-is are not applicable 
to short time average comparisons. 



(Log-Amplitude Covariance) 1 vs. Aperture Diameter . 




& 

g 



p 
0 
P 
0 
6 
P 
• P 

Q 

O 

P 

3 

P 

CD 

PH 

< 

V / 

Q 





V0 


/ — \ 




<NI 


00 


• 


Q 


• 


• 


o 


o 


< 


O 


o 



U 



o 

o 






<4H 

CD 

Pi 



P 

O 

P 

V) 

CD 

•r— I 

"O 

p 

p 

c/) 

p 

0 

> 

CD 

)S 



"P 

O 

♦ i— I 

P 

Ph 

p 

0 

4 -> 

P 



56 



Figure 2 



Scintillation 
Detector Package 




1. Pin 10 Silicon Photodiode 

2. Field Lens 

3. Field Limiting Aperture 

4. Objective Lens 

5. 10 A Narrow-Band Filter 
Centered at 6328 A 

6. 18 Leaf Iris 
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Vertical Position (mm) 



Figure 3 



Pin 10 Photodiode Relative Intensity 
Response (I) vs. Position on Surface. 
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Figure 4 

Scintillation Measurement Electronic Components 
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Figure 5 

Experimental Propagation Path 
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Figure 6 

Laboratory Aperture Averaging Experiment 



im) 




k L___ L 




Sequence 2 



.1 I ! I 1. 






1 



1 




_1 i_ — . 1 - - 1 1 

5 10 IS 20 25 

D(Aperture Diameter, mm) 



61 



Figure 7 

[Log Amplitude Covariance (Using 23 mm Aperture)] 
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Time (POST) 



Figure 8 

1/2 

[Log Amplitude Covariance (Using Aperture D) ] ' 

vs. Time 
22 October 1973 





1600 1700 



Time (P1)ST) 




63 



Figure 9 

1/2 

[Log Amplitude Covariance (Using Aperture D) ] 

vs. Time 
22 October 1973 
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Figure 10 

(23 mm) > vs. Aperture Diameter 
22 October 1973 
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D(Aperture Diameter, mm) 

<C ° a (D) " C^(23 mm) > is the average value of the difference in the 
square roots of the log amplitude variance as sensed through aperture 
diameter D and through a 23 mm aperture. 



[Log Amplitude Covariance (Using 23 mm Aperture)] 

vs. Time 

27 September 1973 



CNJ 




66 



(PDST) 



Figure 12 



[Log Amplitude Covariance (Using Aperture 

vs. Time 
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Figure 13 ]/ 2 

[Log Amplitude Covariance (Using Aperture D) ] 

vs .. T ime 
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<C n ^(D) - C^(23 mm)> is the average value of the difference in the 
square roots of the log amplitude variance as sensed through aperture 
diameter D and through a 23 mm aperture. 



Figure 15 
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D(Aperture Diameter, mm) 

The Square Root of the Normalized Log Amplitude Covariance (D) 

vs. Aperture Diameter 



Figure 16 

1/2 

[Log Amplitude Covariance (Using 23 mm Aperture)] 

vs . 

The Inverse Slope of Regression Lines Plotted in Figure 15. 
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Figure 17 
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Figure 18 
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Figure 19 
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Figure 20 
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Figure 21 
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Figure 22 

1/2 

[Log Amplitude Covariance] vs. Time 
Observed at Two Heights 
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Figure 23 
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Relative Spectral Power [ W ( £ ) / W ( 2 0 H z ) ] 
vs. Frequency 

Observed at two heights above the ocean surface; the upper 
level utilized the variable aperture system, the lower level 
utilized the fixed aperture system (Aperture size, 23 mm). 



78 



Figure 24 
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(6328A) laser beam propagating over a 4.05 kilometer 
horizontal path in the marine boundary layer. The results 
fall in three areas. (1) The covariance of the logarithm of the 
laser beam amplitude increases significantly with decreasing 
detector aperture diameter for apertures less than the first 
Fresnel zone size. This effect appears to be dependent on 
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the magnitude of the refractive index structure constant, 
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C^. (2) The dependence of on (height) for unstable 

meteorological conditions was verified using the laser 
beam in determining C^. (3) Frequency analysis performed 

on scintillation data indicates that the increase of power 
in higher frequencies is more dependent on the wind velocity 
perpendicular to the propagation path than on detector 
aperture diameter for apertures smaller than the first 
Fresnel zone size. 
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